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£-j Abstract 

Wc report the results of our observations of the 12 CO (J= 1-0) and 12 CO (J=3-2) line emission of 74 
major giant molecular clouds (GMCs) within the galactocentric distance of 5.1 kpc in the Local Group 
galaxy M33. The observations have been conducted as part of the Nobeyama Radio Observatory M33 
All-disk survey of Giant Molecular Clouds project (NRO MAGiC). The spatial resolutions are 80pc for 
12 CO ( J= 1-0) and lOOpc for 12 CO (J= 3-2). We detect 12 CO ( J= 3-2) emission of 65 GMCs successfully 
Furthermore, we find that the correlation between the surface density of the star formation rate, which is 
1 derived from a linear combination of Ha and 24/xm emissions, and the 12 CO ( J= 3-2) integrated intensity 

QO . still holds at this scale. This result show that the star-forming activity is closely associated with warm 

P\| ■ and dense gases that are traced with the 12 CO (J— 3-2) line, even in the scale of GMCs. We also find that 

^sO \ the GMCs with a high star-forming activity tend to show a high integrated intensity ratio (i?3-2/i-o)- 

f^*) ■ Moreover, we also observe a mass-dependent trend of i?3_2/i-o f° r the GMCs with a low star-forming 

activity. From these results, we speculate that the i?3_2/i-o values of the GMCs with a low star-forming 
CS| . activity mainly depend on the dense gas fraction and not on the temperature, and therefore, the dense gas 

fraction increases with the mass of GMCs, at least in the GMCs with a low star- forming activity. 
Key words: galaxies: individual(M 33) — ISMxlouds — stars: formation 



$_i ' 1. Introduction Group are preferable, because we can easily construct a 

GMC sample at an uniform distance from us whereas 
The formation processes of massive stars in galaxies are GMCs in the Milky Way are often subject to the distance 
one of the key topics of study in modern astronomy as determination issues. In fact, a wide-area 12 CO ( J= 1-0) 
these massive stars play an essential role in the evolu- survey of the Large Magellanic Cloud (LMC) provides us 
tion of galaxies. Because stars are formed from molecular with rich information about the variation in GMC proper- 
gas, studies on molecular gas in galaxies can provide us ties (Fukui et al. 2008) with the evolution of the associated 
with crucial information on the fundamental physical pro- star clusters (Kawamura et al. 2009). 
cesses of massive star formation. Most of the molecular The disk galaxy M33 is another important target for 
gas is confined to molecular clouds, and virtually all giant such study because while the LMC is an irregular dwarf 
molecular clouds (GMCs) are sites of star formation in galaxy, M33 is one of the nearest galaxies {D = 840 kpc; 
the Milky Way. Thus, the physical properties of the gases Freedman et al. 1991) with a spiral morphology. In ad- 
present in GMCs serve as important information for the dition, its favorable inclination angle (i = 52°; Corbclli & 
investigation of massive star formation. Salucci 2000) enables us to observe the disk face, thereby 

In order to investigate the evolution of GMCs, it is im- making this galaxy unique, 

portant to study GMCs using a statistical approach. For Because of these merits, GMCs in M33 have been statis- 

this purpose, observations of nearby galaxies in the Local tically investigated through observations of the CO lines 
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in several studies using various existing telescopes (e.g., 
for 12 CO(J=1-0): Wilson & Scoville 1990; Engargiola 
et al. 2003; Rosolowsky et al. 2007, and for 12 CO ( J= 2- 
1): Graticr et al. 2012). These studies provide us with 
essential physical properties of GMCs, including molec- 
ular gas masses, sizes, line widths, and virial parame- 
ters. Recent observations using the Herschel PACS and 
SPIRE instruments have investigated the far-infrared and 
sub-millimeter pictures of the cold interstellar medium in 
M33 (HerM33es project: e.g., Kramer et al. 2010). In 
addition to these observation results, understanding the 
formation process of dense molecular gas within GMCs is 
a critical step in understanding the physical connection 
between GMC properties and the evolutionary phases of 
GMCs. Multi-line observations are indispensable tools for 
achieving this understanding. Recent studies have shown 
that the molecular lines with high critical density and/or 
those of high rotational excitation level are closely as- 
sociated with massive star formation in external galax- 
ies. Kohno et al. (1999) found that the HCN emission in 
NGC6951 shows a better spatial correlation with the star- 
forming regions than the 12 CO (J= 1-0) emission; Gao & 
Solomon (2004) reported the existence of tight IR-HCN 
linear correlation in spiral galaxies, which indicates that 
the global star formation rate (SFR) is proportional to the 
mass of dense molecular gas. Mauersbergcr et al. (1999) 
showed the relation between 12 CO ( J= 3-2) emission and 
the SFR at the centers of 28 nearby galaxies. More re- 
cently, the 10-m antenna of Atacama Submillimeter-wave 
Telescope Experiment (ASTE: Ezawa et al. 2004; Kohno 
et al. 2004) was used to conduct a wide-area survey of 
the 12 CO(J=3-2) line toward the inner disk region of 
the nearby starburst galaxy M83 (Muraoka et al. 2007), 
revealing an intimate connection between dense molec- 
ular gas and massive star formation at a spatial scale of 
— 500pc. However, the above mentioned studies are based 
on observations with resolutions greater than a few 100 pc, 
and therefore, observations with higher resolutions are re- 
quired to study the properties of warm and dense gas in 
GMCs. 

Some case studies of dense gas properties in GMCs have 
been conducted for nearby galaxies (Minamidani et al. 
2008; Hirota et al. 2011) including M33; a recent high- 
resolution (~ 100 pc) study of NGC604, the most active 
star- forming region in M33, found an the arclike distribu- 
tion of 12 CO (J=3-2) emission, which is possibly formed 
through compression by stellar winds and/or supernovae 
shocks of previously formed stars (Tosaki et al. 2007). 
Moreover, aperture-synthesis observations of NGC604 by 
Miura et al. (2010) revealed that star-forming efficiency 
decreases as the projected distance from the heart of 
the central cluster increases. Furthermore, they detected 
HCN emission from two massive GMCs at the rim of 
"Ha shells" produced by the cluster. Their result indi- 
cates that the evolutionary stages of GMCs change radi- 
ally in the course of stellar cluster formation. Another 
study of NGC604 using a large velocity gradient model 
for 12 CO(J=1-0), 12 CO(J=3-2), and 13 CO(J=1-0) 
showed that the temperature and density of GMCs vary, 



reflecting their evolutionary stages (Muraoka et al. 2012). 

For conducting statistical surveys of dense molecular 
gas in M33, examination of higher- J CO lines such as 
12 CO (J=3-2) in the emission spectra of GMCs can be 
useful, because these lines are very bright as compared 
with the emission lines of high dipole moment molecules 
such as HCN and HCO+. For instance, two GMCs in 
NGC604 whose HCN emissions were detected through in- 
tcrferometric observations show a HCN-to-CO flux ratio 
(flux unit: Jykms -1 ), S , H CN(i-o)/<5'co(i-o) ! 01 ~l/42 and 
1/56 (Miura et al. 2010). Rosolowsky et al. (2011) re- 
ported no detection of HCN emission from four GMCs in 
M33 at ~ 100 pc resolution, which set very small upper 
limits of the HCN-to-CO intensity ratio (the intensity is 
measured in Kkms -1 ; hereinafter, the same unit is con- 
sidered for intensity), /hcn(i-o)/^co(i-o), from -1/40 
down to 1/280 . Furthermore, for other galaxies, Kuno 
et al. (1995) reported very small /hcn(i-o)/-^co(i-o) ra- 
tios down to —1/50-1/90 in the disk regions of M51, and 
similar ratios were observed for a few 10-pc scale observa- 
tions of GMCs in M31 (Brouillet et al. 2005). In contrast, 
the /co(3-2)/^co(i-o) ratios are usually greater than 1/10 
even in the disk regions of galaxies (e.g., Wielebinski et al. 
1999; Walsh et al. 2002; Muraoka et al. 2007); therefore, 
the 12 CO (J=3-2) line is suitable for whole disk surveys 
of GMCs. 

In this work, we present a study of the warm and 
dense molecular components of GMCs in M33. We in- 
vestigated these components using a combination of a 
high-quality On-the-Fly (OTF) map of 12 CO(J=3-2) 
emission obtained using the ASTE 10-m telescope and 
that of 12 CO(J=1-0) emission obtained by the NRO 
45-m telescope. In this study, we have specifically fo- 
cused only on the physical properties of GMCs; we will 
use a 12 CO(J=3-2) wide-area image and detailed map- 
based studies in a soon-to-be-published paper by Miura 
et al. 2012. This work is part of the Nobeyama Radio 
Observatory M33 All-disk Survey of Giant Molecular 
Clouds (NRO MAGiC) project (Tosaki et al. 2011). 

2. Observations and Data Analysis 

We first selected target regions that cover the majority 
of GMCs in M33, as shown in Fig. 1. Subsequently, we 
conducted 12 CO(J=3-2) observations using the ASTE 
telescope, with a total area coverage of 121 arcmin 2 (~ 
7.3kpc 2 ). The observations were performed from June to 
September 2010 and in November 2011 using the OTF 
observation mode (Sawada et al. 2008). The front-end 
component was a CArtridgc-Typc Sideband-separating re- 
ceiver for the ASTE 345-GHz band (CATS345: Inoue et 
al. 2008). We calibrated the antenna temperature using 
the standard chopper-wheel technique. The main beam 
efficiency was 0.6 ±0.1 at the frequency of the 12 CO ( J= 3- 
2) line. The temperature scale used throughout this paper 
is the main-beam brightness temperature (Tmb ) ■ The typ- 
ical system noise temperature during the observations was 
less than — 200 K in the single sideband (see Miura et al. 
2012 for a more detailed description of the observations). 
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The angular resolution of all images is 25"(~100pc). 
We regridded the images so that their spatial and velocity 
coordinates matched with the 12 CO(J=1-0) data; the 
pixel scale is 7". 5 (^30pc) and the velocity resolution is 
2.5kms- 1 . The rms noise was typically 18-70 mK. 

We used the 12 CO ( J= 1-0) data cube from Tosaki et 
al. (2011). This data was obtained using the NRO 45- 
m telescope and the 25-beam receiver BEARS. A field 
of 30' x 30' (~ 7.3 kpc x 7.3 kpc) was mapped using the 
OTF technique. In that study, the main beam efficiency 
was 0.32 ± 0.02 at the frequency of the 12 CO ( J= 1-0) 
line. The errors in the scaling factors of BEARS were 
smaller than ^20%. The data cube had a velocity width 
of 600 kms -1 , which adequately covers the full velocity 
extent of the emission in the observed area. Its velocity 
resolution was 2.5 kms -1 . The angular resolution and the 
grid spacing were 19". 3 (^80pc) and 7". 5, respectively 
The rms noise was a = 130 mK. 

Wc calculated the SFR using the formulation by 
Calzetti et al. (2007), where the SFR is derived from a 
linear combination of the Ha line and 24-/xm continuum 
emission. We used the Ha image provided by Hoopes & 
Walterbos (2000) and the 24-/xm image from the Spitzer 
Space Telescope data archive. These images were con- 
volved into a 25" beam, regridded to a pixel scale of 7". 5 
to match the 12 CO (J=3-2) data, and finally combined. 
The rms noises in the Ha and 24- fim maps result in a 
Ssfr error of 4.5 x 10~ n M Q yr -1 pc -2 and a sensitivity 
limit of 1.4 x 1O -1O M yr -1 pc -2 . A more detailed de- 
scription of the SFR analyses are presented in Onodera et 
al. (2010). 

We set the criteria for GMC identification from our 
12 CO (J= 1-0) data cube as follows: (i) pixels with the 
12 CO ( J= 1-0) brightness temperature T CO (i-o) should 
be above the noise limit (see below for details) (ii) all pix- 
els in the cube should be conjunct, (iii) the GMCs should 
contain at least three pixels at a velocity channel, and 
(iv) the sizes and velocity widths should be larger than the 
spatial and velocity resolutions, respectively. First, we ap- 
plied the CLUMPFIND algorithm developed by Williams 
et al. (1994) to the 12 CO ( J= 1-0) data cube in order 
to roughly identify the boundaries of the GMCs and ob- 
tain their masses. The CLUMPFIND algorithm requires 
a search increment and the lowest level that is the same as 
the noise limit as an input parameter set. We searched the 
most suitable parameter set to minimize misidentification 
of GMCs and assessed the validity of GMC identification 
by testing many parameter sets and examining the resul- 
tant channel maps of each GMC candidate. In this pro- 
cess, we rejected the parameter sets that resulted in the 
identification of mere aggregations of widespread small 
clumps as single large clump candidates. In such cases, 
the small clumps were not connected with each other in 
either space or velocity domains. We also rejected the pa- 
rameter sets with which most clumps were divided into 
unreasonably small fragments such as those smaller than 
the spatial and velocity resolutions. As a result of this 
parameter search, we adopted the search increment 1.8 a 
and the lowest level 2.7a and roughly identified the GMC 



candidates. Subsequently, we removed the isolated pixels 
that were incorrectly allocated in some of the GMC candi- 
dates. Finally, we rejected the GMC candidates that were 
close to the rim of the data cube as well as the ones whose 
sizes and velocity widths were smaller than the spatial and 
velocity resolutions. 

We successfully determined 74 GMCs that were also 
contained in our target fields of 12 CO (J= 3-2). Their 
galactocentric distance ranges from 0.2 to 5.1 kpc. The H2 
gas mass of each GMC is calculated from the 12 CO ( J= 1- 
0) data assuming the 12 CO-to-H2 conversion factor Xqo = 
2 x 10 20 cm- 2 (Kkms -1 )- 1 (Rosolowsky et al. 2003). The 
coordinates and velocity at the Tco(i-o) peak position 
and the H2 gas mass of these GMCs are presented in Table 
1. 

3. Results 

3.1. 12 CO(J=1-0) and u CO(J=3-2) line spectra 

Fig. 2 shows the 12 CO(J=3-2) and 12 CO(J=1-0) 
spectra at the Tcori-o) peak positions of each GMC. 
The 12 CO(J=1-0) images were convolved into a 25" 
(~ 100 pc) beam and regridded to a pixel scale of 7". 5 
in order to match the 12 CO (J=3-2) data. Our criterion 
for detection in 12 CO(J=3-2) is that the spectra show 
emission above the 3-sigma level. The rms noise level for 
each spectrum is shown in Fig. 2. 

Among the 74 GMCs, nine arc not detected for 
12 CO(J=3-2) emission (Nos. 25, 51, 52, 56, 62, 63, 64, 
72, 73). Although some of these spectra show multiple 
peaks, the velocity components of the sub-peaks are ex- 
cluded from the GMCs in the process of identification. 
Thus, each GMC corresponds only to the main peak in 
each spectrum. We fitted the main peaks in these spectra 
with a single Gaussian to obtain the integrated intensi- 
ties of CO lines at the peak positions of the GMCs. The 
integrated intensities are listed in Table 1. 

3.2. 12 CO (J =3-2) to 12 CO (J =1-0) integrated inten- 
sity ratio i? 3 _2/i-o 

A tight correlation exists between the 12 CO ( J= 3-2) in- 
tegrated intensity, ico(3-2)> an d the 12 CO (J= 1-0) inte- 
grated intensity, /co(i-o)i m nearby galaxies at a ^500pc 
resolution (Muraoka et al. 2007). Fig. 3 shows the plot 
of /co(3-2) versus ico(i-o) 01 the GMCs in M33. We 
found that the correlation between ico(3-2) an d ^co(i-o) 
still holds at the GMC scale with the correlation coeffi- 
cient of 0.65, whereas i? 3 _ 2 /i-o = ^co(3-2)/^co(i-o) is 
spread over a wide range from less than 0.1 to 0.74. When 
weighted by the measurement uncertainty, the i?3_2/i-o 
mean is 0.26. The GMCs that are not detected for 
12 CO (J=3-2) emission are not used for averaging. This 
value of i? 3 _2/i_o is fairly low as compared with the disk- 
averaged 12 CO ( J= 2-l)/ 12 CO (J= 1-0) of 0.73 discussed 
by Gratier et al. (2010). This indicates that the tempera- 
ture and local volume density at the GMC center is gener- 
ally not very high at the resolution of ~ 100 pc. Although 
the GMCs we detected are predominantly located in spi- 
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ral arms, the i?3_ 2 /i_o value is slightly smaller than that 
of the arm region of the Milky Way galaxy (~0.4, Oka et 
al. 2007). A CMC with a remarkably high i? 3 _ 2 /i-o value 
of 0.74 (No. 13) is located at the center of the most active 
star forming region, NGC604. Such elevated R3-2/1-0 
values are also reported in the most active HII regions in 
LMC, 30 Dor, and N159 (i? 3 - 2 /i-o -0.7, Minamidani et 
al. 2008), and in larger-scale measurements of the nearby 
starburst galaxy M83 (-0.7, Muraoka et al. 2007), lumi- 
nous IR Galaxies (LIRGs, —0.6, Yao et al. 2003), and 
starbursting compact dwarfs (—0.7, Israel 2005). 

4. Discussion 

In this section, we focus on the physical properties of 
our sample GMCs and their relationship with massive star 
formation. In particular, we discuss the relationship be- 
tween the denser section of molecular clouds and massive 
star formation on a scale of — 100 pc, as well as the possi- 
bility of the existence of a high dense-gas fraction within 
massive GMCs. Here, we used the 12 CO (J~ 1-0) and 
the SFR data that were convolved into a 25"(— lOOpc) 
beam so that their resolution was the same as that of the 
12 CO(J=3-2) data. 

4-1. Relationship between SFR and CO intensities 

Fig. 4 shows the plot of Ssfr against /co(i-o) an( l 
Jco(3-2) at the Ico(i-o) peak position in each GMC. It 
should be noted that Ssfr is integrated over all the veloc- 
ity components in the line of sight, whereas the /co(i-o) 
and ^co(3-2) values correspond only to the velocity ex- 
tent of the GMC. Because the sub-components in the CO 
spectra of most GMCs are smaller than the main compo- 
nent, this effect on Ssfr is estimated to be a factor of less 
than 2. 

This figure shows a clear correlation between /co(3-2) 
and Ssfr, with the correlation coefficient being 0.68. In 
contrast, the correlation between Ico(i-o) an d Ssfr is 
not obvious, with the correlation coefficient being 0.22. 
This result is consistent with that obtained by Yao et al. 
(2003); their study showed that the far- infrared (FIR) lu- 
minosity and the 12 CO (J= 3-2) luminosity were better 
correlated than the FIR luminosity and the 12 CO ( J= 1- 
0) luminosity in infrared-luminous galaxies with resolu- 
tions of — 1 — 14 kpc. Our result is also consistent with 
that of Komugi et al. (2007); they presented an identical 
result for the nuclei of dozens of nearby galaxies, with a 
typical resolution of — 3 kpc, and obtained a correlation 
in the ico(3-2) range of —10 — 10 3 Kkms~ 1 . We found 
that this correlation is valid down to the ^co(3-2) range 
of —1 — lOKkms -1 , as shown in Fig. 4, even in the GMC 
scale of — 100 pc. 

This correlation can be attributed to the higher den- 
sity and/or temperature of the gas that is traced with the 
12 CO(J=3-2) line. That is, because the 12 CO(J=3-2) 
line traces denser gas (with a critical density of n cr j t — 
10 4 cm -3 ) than that traced by the 12 CO(J=1-0) line, 
with n cr ;t — 10cm -3 , a good correlation is expected, as 
seen between the HCN intensity and the SFR, even up 



to the dense core scale (Wu et al. 2005). Furthermore, 
star-forming activities heat the surrounding gas, and con- 
sequently, the CO molecules in the gas are excited to show 
stronger J~ 3-2 line emission. Because the CO J = 3 level 
is 33 K above the ground state and the CO J = 1 level is 
5.5 K, we can also expect that the correlation between the 
SFR and the 12 CO(J=3-2) line emission is better than 
that between the SFR and 12 CO ( J= 1-0) line emission. 

As we have shown in the previous subsection, the in- 
tegrated intensities of 12 CO (J= 1-0) and 12 CO (J=3-2) 
are well correlated with each other, at the peak positions 
of Tco(i-o) m each GMC, at a resolution of — lOOpc. 
Further, considering that Ssfr and ^co(3-2) are better 
correlated than Ssfr and /co(i-o)i we expect a radial 
gradient of emissions of the CO lines and the SFR trac- 
ers in the GMCs. From the center of the GMCs to that 
of the HII regions, we observe peaks of 12 CO (J= 1-0), 
12 CO(J=3-2), and Ha+24^m. The existence of the 
radial gradient can be inspected by directly comparing 
the maps of 12 CO (J= 1-0), 12 CO (J= 3-2), and the SFR 
tracers. The results of this comparison will be presented 
in our forthcoming paper. 

4-2. Relationship o/i? 3 _ 2 /i-o with GMC mass and SFR 

As mentioned above, the emission of the 12 CO (J=3 
2) line is more closely related to the SFR than that of 
the 12 CO(J=1-0) line is. Moreover, it is probable that 
i? 3 _2/i-o depends on the properties of the GMCs, such 
as their mass. For example, in our galaxy, clouds with 
masses greater than 10 3 Mq have denser and warmer cores 
than those with masses less than 10 3 Mq (Vcrschuur & 
Kellermann 1988). 

Fig. 5 shows a plot of the H 2 gas mass of the GMCs 
against i?3_ 2 /i_o towards the Tco(i-o) peak position of 
each GMC. In this figure, the GMCs are denoted by cir- 
cles, whose area is proportional to the Ssfr at the same 
position. The red circles indicate the GMCs with a high 
star-forming activity (Ssfr >1x 10 -8 Mq pc -2 ), and the 
blue ones indicate those with a low star-forming activ- 
ity (S S fr < 1 x 10 -8 Mq pc -2 ). Here, in order to check 
whether R 3 _ 2 /i-q shows mass-dependency, we divide the 
sample in half at Mh 2 = 8.4 x 1O 4 M . The crosses indicate 
the averages of i?3_ 2 /i_o m each mass bin; these averages 
are taken for GMCs with high and low star-forming ac- 
tivities. The averages of -R3-2/1-0 are 0.31 for the 14 
less-massive and star-forming GMCs, 0.37 for the 23 mas- 
sive and star-forming GMCs, and 0.24 for the 14 mas- 
sive and low star-forming GMCs. It should be noted that 
R3-2/1-0 = 0.14 for the 23 less star-forming, less-massive 
GMCs is a value averaged over upper-limit values. 

We found that the GMCs with high star-forming ac- 
tivity tend to show a high i? 3 _ 2 / 1 _ value. Moreover, 
the mass-dependent trend of i?3_ 2 /i_o m the GMCs with 
a low star-forming activity is apparent. We found that 
- 24% of the GMCs with M H2 < 8.4 x 10 4 M show 
i?3-2/i-o < 0.1, whereas none of the GMCs with Mh 2 > 
8 . 4 x 1 4 Mq have such a small value of i? 3 _ 2 / 1 _ ■ This dif- 
ference between the massive and less-massive GMCs is the 
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reason for the appearance of a blank in the lower-right por- 
tion of this plot. It also appears in the averaged i?3-2/i-o : 
the difference between the mean upper- limit value of 0.14 
for less star-forming, less-massive GMCs and the mean 
value of 0.24 for less star-forming, massive GMCs is sig- 
nificant: it is larger than the measurement uncertainties. 
Because the Esfr in the less star-forming GMCs does not 
greatly depend on GMC mass, the mass-dependent trend 
of -R3-2/1-0 is attributed to factors other than Ssfr- On 
the other hand, -R3-2/1-0 m the GMCs with a high star- 
forming activity does not significantly differ between the 
massive and less-massive GMCs: the difference between 
the average values of i?3_2/i-0j i-e-, 0.37 and 0.31, is com- 
parable with the measurement uncertainties. 

On the basis of these results, we speculate that the val- 
ues of -R3-2/1-0 for the GMCs with a low star- forming 
activity mainly depend on the dense gas fraction and 
not on the temperature, because they have a relatively 
small amount of energy input from the heating source. 
Considering the mass-dependent trend of the R3-2/1—0 
values of the GMCs with a low star-forming activity, we 
can state that the dense gas fraction increases with the 
GMC mass, at least in the case of GMCs with a low 
star-forming activity. Because the SFR and dense gas 
mass show a linear correlation (Wu et al. 2005), the mass- 
dependency of the dense gas fraction is consistent with the 
results obtained by Kennicutt et al. (2007), who showed 
that the star-formation efficiency increases (or the star for- 
mation timescale reduces) with an increase in the cloud 
mass. On the other hand, i? 3 _ 2 /i_o of the GMCs with 
star-forming activity does not significantly differ between 
the massive and less-massive GMCs, because it depends 
mainly on the temperature. 

5. Summary and Conclusions 

We identified 74 GMCs from our 12 CO ( J= 1-0) data 
of M33. Among these, 12 CO(J=3-2) emission of 65 
GMCs are detected. We compared the 12 CO (J= 1-0) and 
12 CO (J=3-2) line emissions of the GMCs to investigate 
the variation in their physical properties and evolutionary 
stages. 

1. The 12 CO(J=1-0) and 12 CO(J=3-2) integrated 
intensity ratio ^3-2/1—0 is spread over a wide range 
from less than 0.1 to 0.74, having a weighted mean 
of i?3_2/i-o = 0.26. This weighted mean is slightly 
smaller than that of the quiescent disk region of the 
Milky Way. 

2. The surface density of SFR and the 12 CO {J =2,-2) 
integrated intensity at the 700(1-0) peak position 
of GMCs show a good correlation, which is consis- 
tent with the findings of previous studies that used 
data of lower resolutions (^500pc). We have shown 
that the correlation is valid down to /co( J = 3-2) ~1- 
lOKkms" 1 , even in the scale of GMCs. 

3. We have investigated the relationship among the 
GMC mass, i? 3 _ 2 /i_o, and the SFR. We have found 
that the GMCs with a high star-forming activity 



tend to show a high -R3-2/1-0 value. We have also 
observed the mass-dependent trend of -R3-2/1-0 in 
the GMCs with a low star-forming activity. On 
the basis of these results, we speculate that the 
-R3-2/1-0 values of the GMCs with a low star- 
forming activity mainly depend on the dense gas 
fraction and not on the temperature. This implies 
that the dense gas fraction increases with the mass 
of GMCs, at least in the GMCs with a low star- 
forming activity. 

We would like to thank the anonymous referee for his/her 
comments, which have helped enhance this paper. The 
Nobcyama Radio Observatory is a branch of the National 
Astronomical Observatory of Japan, National Institutes 
of Natural Sciences. 
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Table 1. Observed properties of GMCs 



No. 




RA 


Decl 


Vrad 


/co(j = i- 


-0) 


^CO(J=3-2) 


R31 


GMC mass Af H , 


10 a x S SFR 






(J2000) 


(J2000) 


(kms- 1 ) 


(K kms" 


') 


(Kkms- 1 ) 




(10 4 M Q ) 


(M Q yr- 1 pc- 2 ) 


(1) 




(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


1 


1 


34 09.65 


30 49 07.5 


-248.8 


21.53 




5.66 


0.26 ± 0.02 


66.3 


3.2 


2 


1 


34 13.66 


30 33 45.0 


-157.2 


12.33 




6.41 


0.52 ± 0.07 


20.9 


36.9 


3 


1 


34 06.73 


30 47 45.0 


-256.7 


11.34 




3.14 


0.28 ± 0.03 


24.1 


22.5 


4 


1 


33 35.90 


30 39 30.0 


-167.9 


9.36 




4.36 


0.47 ± 0.05 


38.5 


27.9 


5 


1 


33 49.85 


30 33 52.5 


-133.7 


7.46 




2.64 


0.35 ± 0.04 


16.2 


13.3 


6 


1 


34 10.77 


30 36 15.0 


-159.3 


7.80 




3.41 


0.44 ± 0.03 


13.8 


20.8 


7 


1 


33 59.75 


30 49 15.0 


-252.0 


10.31 




3.75 


0.36 ± 0.03 


17.3 


6.0 


8 


1 


33 35.91 


30 36 30.0 
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9.81 




3.30 


0.34 ± 0.04 


26.3 


20.6 


9 


1 


33 41.14 


30 39 15.0 


-165.8 


14.89 




4.68 


0.31 ± 0.03 


35.6 


8.9 


10 


1 


34 34.66 


30 46 15.0 


-221.2 
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2.84 


0.24 ± 0.03 


14.2 


17.4 


11 


1 


33 58.58 


30 48 52.5 


-245.9 


11.98 




4.18 


0.35 ± 0.04 


33.3 


11.5 


12 


1 


34 08.45 


30 39 15.0 


-195.0 


10.75 




3.24 


0.30 ± 0.03 


21.4 


18.1 


13 


1 


34 33.50 
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-241.7 
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0.74 ± 0.07 


40.9 


263.9 
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33 29.54 
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27.2 
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2.47 


0.21 ± 0.02 


17.8 


2.2 


17 


1 


33 44.63 


30 36 00.0 


-139.3 


8.17 




3.85 


0.47 ± 0.06 


18.2 


17.3 


18 


1 


34 02.64 
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22 


1 


33 52.76 


30 39 15.0 
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8.51 




5.01 
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11.4 


48.3 


23 


1 


33 47.53 


30 32 37.5 
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9.97 




2.69 
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15.9 


10.3 


24 


1 


33 52.76 
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-154.8 


7.11 




2.33 
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6.4 


25 


1 


34 14.25 
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11.8 


2G 
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33 38.81 


30 41 00.0 


-186.9 
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1 


34 02.06 


30 38 37.5 
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Table 1. (Continued.) 
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RA 
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1 
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0.79 
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1 
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GO 


1 
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16.8 


Gl 


1 
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1.72 


0.27 ± 0.06 


2.3 


23.6 


G2 


1 


34 04.38 


30 39 22.5 
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3.70 


<0.20 


<0.05 


2.2 


3.9 


G3 


1 


34 10.76 


30 35 07.5 


-155.5 


2.68 


< 1.34 


<0.50 


2.5 


3.5 


G4 


1 
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-242.6 


5.28 
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1.6 


4.6 


G5 


1 


34 13.10 
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-197.0 


4.21 
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2.3 


2.7 


GG 


1 


34 07.31 


30 47 52.5 


-256.2 


10.21 


2.93 


0.29 ± 0.05 


2.8 


11.6 


G7 


1 


33 52.76 
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-209.0 


4.05 


0.36 
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2.3 
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G8 


1 
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0.72 
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4.0 


3.3 


69 


1 
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4.13 


0.52 
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2.1 


5.0 


70 


1 


34 14.88 


30 46 22.5 


-249.6 


3.03 


0.40 


0.13 ± 0.04 


2.8 


3.3 


71 


1 


33 48.69 


30 39 45.0 


-174.5 


10.85 


3.26 


0.30 ± 0.05 


2.3 


29.0 


72 


1 


33 30.12 


30 32 45.0 


-124.8 


2.46 


<0.40 


<0.16 


1.7 


3.0 


73 


1 


33 21.40 


30 34 00.0 


-131.7 


5.08 


<0.70 


<0.14 


2.4 


0.7 


74 


1 


33 45.79 


30 33 22.5 


-126.2 


3.82 


1.76 


0.46 ± 0.13 


1.7 


10.8 


NOTES. 


—Col. (1): 


Running Number. Col. (2), (3) and (4): The peak position of T co{1 „ 0) in GMCs. Col.(5):The "CO ( J = 


1-0) integrated 



intensity on a Tmb scale derived by using a single Gaussian fitting for a spectrum at the Tqo(i-o) peak position in GMCs. The data were 
convolved into a 25"(~100pc) beam to match the resolution of the 12 CO(J=3-2) data. Col. (6) Same as col. (5), but for 12 CO(J=3-2) 
line. Col. (7): Ratio of / C 0(3-2) to -f C O(i-0)- Col.(8): The H 2 gas mass of GMCs. Col.(9): * For GMCs whose 12 CO(J=3-2) emissions 
are not detected, ^cO(J = 3— 2) ' s shown as the upper limit, which is estimated as three times the rms of ^cO(J = 3-2) that is calculated 
from the rms noise of Tqo(j = 3_2) an d the width of 12 CO (J= 1-0) line. 
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Fig. 1. Our target regions of the 12 CO (J=3— 2) line drawn on the 12 CO (J= 1-0) integrated intensity gray-scale map of M33 by 
Tosaki et al. (2011). 
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Fig. 2. 12 CO(J=1-0) (black), and 12 CO(J=3-2) (red) spectra and fitted Gaussians of each GMC on a T MB scale. The 
12 CO(J=1-0) images were convolved into a 25" (~ 100 pc) beam and regridded to a pixel scale of 7". 5 in order to match the 
1 CO(J=3-2) data. The running number of the GMCs is shown in the upper left and the 1-sigma level of the 12 CO (J=3-2) 
spectra in a unit of K in the upper right of each panel. Gaussians are not fitted for the 12 CO (J= 3-2) spectra in which significant 
emission is not detected. 
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Fig. 3. Comparison of the 12 CO (J=3-2) and 12 CO (J= 1-0) integrated intensities (on a Tmb scale) at the Tqq(i— o) peak position 
in each CMC. The 12 CO(J=1-0) data were convolved into a 25"(~100pc) beam to match the resolution of the 12 CO(J=3-2) 
data. The line indicates the weighted-mean value of i? 3 _2/i_o = 0.26. The GMCs that are not detected for 12 CO ( J= 3—2) emission 
are not used. 
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Fig. 4. The surface density of SFR (Esfr) versus icO(3— 2) and ^CO(l— 0) at * ne ^CO(l— 0) peak position in each GMC. The 
12 CO ( J= 1-0) and the SFR data were convolved into a 25"(~100pc) beam to match the resolution of the 12 CO ( J= 3-2) data. The 
dashed lines indicate the sensitivity limits for each CO line. As mentioned in Section 2, the error and the sensitivity limit of Esfr. 
are 4.5 X 10~ n M Q yr^pc" 2 and 1.4 X 10~ 10 M Q yr^pc" 2 , respectively. The GMCs that arc not detected for 12 CO(J=3-2) 
emission are not used. 
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Fig. 5. GMC mass versus i?3_2/i— 0- The area of the circles is set proportionally to Esfr at the Tqo(i— o) peak posi- 
tion of each GMC. The red circles represent GMCs with Ssfr > lx 10~ 8 Mq pc -2 , and the blue ones indicate those with 
EgpR < 1 X 10~ 8 Mq pc~ 2 . The dotted line indicates the criterion of Afjj2 = 8.4 X 10 4 Mq. The black crosses and the num- 
bers next to them stand for the averaged i?3_2/i_o values in each mass-bin, which are taken for GMCs with higher and lower 
star-forming activity, respectively. The numbers in the bracket are the numbers of the GMCs used for averaging. The R3—2/1—0 
mean for the less star-forming, less-massive GMCs is an upper- limit value. The 12 CO(J=1-0) and the SFR data were convolved 
into a 25"(~100pc) beam to match the resolution of the 12 CO(J=3— 2) data. For GMCs whose 12 CO(J=3— 2) emissions are not 
detected, R3—2/1— ls shown as the upper limit, with an arrow whose length is the rms value of -R3_2/i— 0- 



